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Transverse modes and instabilities of a bunched beam 
with space charge and resistive wall impedance 



(N 

o 

(N 

C 

c3 



V. Balbekov 
Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, Illinois 6051^\ 
(Dated: January 25, 2012) 

Transverse instability of a bunched beam is investigated with synchrotron oscillations, space 
charge, and resistive wall wakefield taken into account. Boxcar model is used for all-round analysis, 
and Gaussian distribution is invoked for details. The beam spectrum, instability growth rate and 
effects of chromaticity are studied in a wide range of parameters, both with head-tail and collective 
bunch interactions included. Effects of the internal bunch oscillations on the of collective instabilities 
is investigated thoroughly. Landau damping caused by the space charge tune spread is discussed, 
and the instability thresholds of different modes of Gaussian bunch are estimated. 

PACS numbers: 29.27.Bd 
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I. INTRODUCTION 

Transverse instability of a bunch in a ring accelera- 
tor has been considered independently by Pellegrini [1| 
and Sands with synchrotron oscillations and some in- 
ternal degrees of freedom of the bunch taken into ac- 
count ("head-tail instability"). Later Sacherer have in- 
vestigated the effect in depth including dependence of 
the bunch eigenmodes on amplitude of synchrotron os- 
cillations ( "radial modes" Q). 

Space charge field of the bunch was not taken into ac- 
count in these pioneer works. Its role was studied first in 
Ref. [H where it has been shown that the space charge 
tune spread results in Landau damping which suppresses 
many of the head-tail modes, much like other sources of 
the incoherent tune spread. 

However, the last result has been obtained with an as- 
sumption that the space charge tune shift is significantly 
less than the synchrotron tune. A closer examination of 
the problem in Ref. [|| has led to the conclusion that 
almost all head-tail modes are prone to Landau damp- 
ing till then the space charge tune shift is about twice 
less then the synchrotron tune. The damping vanishes 
when the shift becomes more, lower eigenmodes being 
free from the decay first. The lowest (rigid) mode is the 
only universal exception from the rule being potentially 
unstable with any space charge. Sometimes 1-2 next 
modes can demonstrate similar behavior, in dependence 
on the bunch shape. 

However, no wakefield was included in the analysis of 
Ref. [5[ in fact, so the results might be interpreted only 
as conditions which permit the bunch instability but do 
not determine its characteristics completely. Effects of 
short wakes were actually studied in Ref [6fl where the 
instability growth rate has been found at different bunch 
parameters like its length, chromaticity, etc. Transverse 
modes coupling instability was considered in the work as 



well, and it has been shown that this more rigorous effect 
can be caused only by positive wake (for comparison, 
resistive wall creates a negative wake) . 

A single bunch was investigated in the mentioned arti- 
cles because only short-range forces have been accepted 
in all the cases. More general theory is developed in 
this paper where a bunched beam with arbitrary num- 
ber of the bunches is examined taking into account the 
space charge, intra-bunch and bunch-to-bunch interac- 
tions. Both kinds of the interactions affect the beam 
eigenmodes and the instability growth rate, although one 
of them can dominate is specific situation. Just from this 
standpoint the common used terms like "collective modes 
instability" or "head-tail instability" are treated in the 
paper. 

The presentation is focused first on the resistive wall in- 
stability Q-Jl] but the results can be rather easy adapted 
for other known wakefields. The boxcar model is inten- 
sively used to get a general outlook of the problem in a 
wide range of parameters. More realistic Gaussian bunch 
is closely examined in a limiting case of low synchrotron 
frequency, although the opposite case of high frequency is 
also invoked to estimate thresholds of Landau damping. 

Both betatron and synchrotron bare oscillations are 
implied to be linear in this paper. The nonlinear effects 
unquestionably require a special investigation, especially 
as an additional factors of Landau damping. 



II. BUNCHED BEAM GENERAL EQUATIONS 

With space charge and wakefield taken into account, 
equation of coherent betatron oscillations of a bunched 
beam in the rest frame can be written in the form |Q|: 
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where 9 is longitudinal coordinate (azimuth), to is fre- 
quency of the coherent oscillations, Qo an d are cen- 
tral betatron tune and revolution frequency (in the labo- 
ratory frame), AQ(9) is the space charge tune shift aver- 
aged over transverse directions, Q s and <p are synchrotron 
tune and phase, C — ~£,/ 7 l is normalized chromaticity, 
and q{9) is normalized wakefield which specific form will 
be defined later. The bunch shape in the longitudinal 
phase space (9, u) is described by a distribution function 
F(9, u) with corresponding linear density: 
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(normalization of the functions will be specified later). 
The function Y(9,u) is an average transverse displace- 
ment of particles located in the point (9, u) of the longi- 
tudinal space, multiplied by the factor exp [i (Qq + 
Additional averaging of this function over momentum is 
denoted as Y(9) being defined by the expression: 



p(6)Y(t,d) 



F(9,u)Y(t,9,u) du 
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It is more convenient to represent the right hand part 
of Eq. (1) as a sum over all the bunches preceding the 
examined one (including all preceding turns): 
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where v — uj/CIq — Q , A = £ — v ~ £, M is num- 
ber of bunches in the beam, p m {9') — p(2nm/M + 9'), 
and Y m (9') = Y(2nm/M + 9'). Any integral in this ex- 
pression is taken over one of the bunches which are not 
presumed yet to be identical (in particular, some of them 
may be empty). However, periodicity conditions must be 
satisfied in any case: 



Pm+id{Q) = Prn(O), 

Y m+M {9) = Y m {9) exp [-2tu(Q 
A. Symmetric beam 
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In subsequent sections we will consider symmet- 
ric beams composed of identical equidistant bunches: 
Pn{9) = p{9). It would be checked than that all the 
solutions of Eq. (4) fall into M groups which are known 
as collective modes: 
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where k — 1,2,...,M. There are the head-tail modes 
inside the groups, each satisfying the equation: 
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where 29$ is the bunch length, and k = (k—ui/flo)/M ~ 
(k — Qq)/M. Note that the head-tail modes are not au- 
tonomous formations but each of them more or less de- 
pends on the collective mode which it falls. 

The variable t = 9/9$ will be used further as a new 
longitudinal coordinate having a range [—1,1], and the 
normalization condition will be imposed 



J p(r)dr=l. 



III. RESISTIVE WAKE 
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Resistive wall instability was predicted first it Ref. [3]- 
Q. Corresponding wakefield function q(9) is negative, 
and can be presented in the convenient form: 
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where Nb is number of particles per bunch, ro = 
e 2 /mc 2 is classic radius of the particle, R is the machine 
radius, a is specific conductivity of the beam pipe, b y is 
its semi-height, and a is the pipe form- factor (a = 1 for 
a round pipe). Then one can rewrite Eq. (7) in the form: 
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where T = 2tt/M9q is the bunch spacing in terms of 
the variable r, and the notation is used for a shortness: 
y(r) = Y(t)p(t) exp(— iXOqt). Square root in the last 
integral not much depends on the addition (r' — r), so 
the expansion into Taylor series can be applied resulting 
in: 

!>r AQ (Y - Y) ~ q o exp(iA0 o r) 
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FIG. 1: Functions w^c|Vi(fi;) and Via(ft) (resistive wall wake). 



FIG. 2: Eigentunes of the boxcar bunch 



where B = MOq/tt is bunch factor, and 
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Both the functions are periodical, one period of V2(k) be- 
ing plotted in Fig. 1. Because the function Vi(k) has 
the singularities, its plot is multiplied by aperiodic factor 
■\/\k\ for a convenience. 

It goes without saying that form of these plots depends 
on the wakefield. However, in most cases it does not 
change a general structure of Eq. (11) as well as many 
subsequent inferences. 



IV. BOXCAR MODEL 

The boxcar model is described by the expressions 
1 , , 1 
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where u is normalized transverse momentum conjugated 
with the longitudinal coordinate r: u 2 = A 2 — r 2 , A = 
T max- A virtue of this model is that all solutions of 
Eq. (11) are exactly known at go = with any param- 
eters Q s and AQ, as it was shown first in Ref. [3J and 
developed in detail in Ref. [5]. This circumstance gives 
a great possibility to overlook the wakefield produced ef- 
fects. Required information is shortly reminded below. 

All the solutions of Eq. (11) with q = are derivable 
from Legendre polynomials Y(t) = P„(r), n = 0, 1, 2, ... 
At any n, there are n + 1 different eigenfunctions 



Yn,m{ r i u ) with eigentunes v n<m where m = n,n — 
2, 2 — n, — n. They satisfy the orthogonality condi- 
tions: 



(15) 



All the eigentunes are real numbers. Some of them are 
plotted in Fig. 2 where the space charge tune shift is 
used for a scaling. At fj, = Q S /AQ — 0, the eigen- 
tunes f n<n take start from the point v — that is 
lj/Qq = Qo (the bare tune). Other eigentunes v n ,m<n 
start from the point v = —AQ which corresponds the 
incoherent tune with the space charge included: w/^o = 
Qo — AQ. At [i <C 1, all eigentunes acquire the ad- 
ditions 8v oc Q 2 /AQ. In this limiting case, the eigen- 
functions have about linear polarization: Y n<n depends 
mostly on longitudinal coordinate r at any n, and all 
other functions Y n , m<n depend mostly on momentum 
u (many examples are given in Ref. 0). However, 
these solutions merge together at /i> 1, forming multi- 
poles Y n , m {A) exp(iTO0) with different dependence from 
synchrotron amplitude. These radial modes Y n>m {A) are 
born by Legendre polynomials of powers n = |m|, |m| + 2, 
etc. Note that the functions Y n ^ n (A) are the lowest (min- 
imally oscillating) radial modes with given n = to. Their 
eigentunes are n(n + l)/2 x Q 2 /AQ at small [i and 
~ nQ s at large one. The functions Y n ^ m<n are treated 
usually as higher radial modes. 



A. Low wake 

The study of the boxcar model is continued in this 
subsection with an assumption that the wake is small 
enough to apply the perturbation methods (applicability 
of this approximation will be discussed in Sec. V). Then 
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FIG. 3: Functions A„ jm (/i) of the boxcar bunch. Numbers 
(n,m) are given near each line. 



FIG. 4: Functions Re / n (#oC) °f the boxcar bunch. The right- 
hand parts of these even functions are shown. 



A = £ in Eq. (11), and the additions to the eigentunes 
can be presented in the form: 
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Generally, these relations are applicable to a bunch of 
any form with low wakefield. However, the boxcar model 
is the only known case at present which allows to inves- 
tigate the problem in depth because its eigenfunctions 
Y m , n (T,u) andy„(r) = P„(r)p(r) exp(-i9 (r) are really 
known with any \i. The results are presented graphically 
in Figs. 3-7 and commented below. 

The coefficients A„ j7n depend only on the ratio \x = 
Q S /AQ, and does not depend on the bunch length or 
chromaticity (Fig. 3). They describe a general effect of 
synchrotron oscillations and space charge on transverse 
coherent motion of the bunches, including possible insta- 
bility growth rate. It is seen that the space charge tune 
shift enhances influence of the wakefield on the lowest ra- 
dial modes Y nn but depresses its influence on the higher 
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FIG. 5: Functions Im /n(#oC) of the boxcar bunch. The right- 
hand parts of these odd functions are shown. 



modes Y n m<n . It is very understandable result because 
higher radial modes are polarized mostly in ^-direction at 
small /i, so that the global bunch displacement at given 
azimuth and, correspondingly, excited wakefield should 
be relatively small in this limiting case. 

In contrast with this, the part of Eq. (16) in square 
brackets describes effects of the bunch length and chro- 
maticity on different n-modes. There are three terms 
here which are concerned with different physical effects. 

Interaction of particles inside the bunch is described 
by the coefficients / n (#oC) some of which are plotted in 
Figs. 4-5. With non-zero chromaticity, this part is capa- 
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FIG. 6: Functions <7 n (#oC) °f the boxcar bunch. The right- 
hand parts of these even functions are shown. 
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FIG. 7: Functions h n (9oQ of the boxcar bunch. The right- 
hand parts of these even functions are shown. 



ble to cause the head-tail instability of different modes 
which basic properties were predicted in earliest works 
[l[-0. It is a single-bunch effect which is proportional 
to the bunch population and does not depend on number 
of the bunches. 

Second term in Eq. (16) describes the main effect of 
the collective interaction of the bunches (factor Vi(rt) in 
Fig. 1) including its dependence on chromaticity (factor 
9n(@o£))' As a rule, the term gives a maximal contribu- 
tion to the tune shift, especially if the beam consists of 
many bunches. Indeed, with M S> 1 one can get k - 
y/(k-Qo)/M < 1 that is V^k) ~ (1 + k/MWR 
for the most unstable modes (k > Qq). Then the total 
contribution of this term to the tune is: 

q MA 
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A^n,m ^ ,, 1+ I T, 7T7 ( 21 ) 

With Eq. (9) for q$ used, it gives the expression: 
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where N — MNj, is the total beam intensity, uj = 
fio(& — Qo) is the coherent frequency in the laboratory 
frame, and 8 (u>) is corresponding skin depth. First part 
of the formula coincides with well known expression for 
resistive wall instability of a coasting beam Q-fl]. The 
factors A„ !m describe an impact of the bunching upon 
different head-tail modes, including their dependence on 
synchrotron frequency and space charge tune shift. Note 
that the coefficient of the most important rigid mode 
Ao,o — 1 independently on the mentioned parameters, 
that is the bunching does not affect this mode. 

The last term g n (6o() in Eq. (22) describes an influence 
of chromaticity. In this regard, it is pertinent to dwell 



upon the different character of the chromatic effects in 
coasting and bunched beams. 

In first case, chromaticity leads to a spread of incoher- 
ent betatron frequencies which phenomenon can cause 
Landau damping resulting in total suppressions (preven- 
tion) of instability. 

In contrast with this, average momenta of all the par- 
ticles are equalized in the bunch through synchrotron os- 
cillations. In such conditions, chromaticity does not pro- 
vide a systematic tune spread and cannot bring about 
Landau damping at once. There is no questions that 
additional slip of betatron phases of particles with re- 
spect to the coherent field phase affects the interaction 
and can change the instability growth rate. As it follows 
from Fig. 6, switch of the instability peaks to the higher 
internal modes is the most descriptive result of this. 

However, it is necessary to take into account also other 
then chromaticity factors which can result in Landau 
damping. In particular, space charge tune spread itself 
can produce this effect in bunched beams. As it is shown 
in Ref. JU-JU, the higher internal modes are more sensi- 
tive to this kind of Landau damping. Therefore, one can 
expect that the above mentioned shift of the instability 
peaks is capable to suppress all internal bunch modes ex- 
cept the rigid one which is not prone to this kind of the 
damping [5J. However, this problem cannot be solved in 
frames of the boxcar model which ignores this part of 
the tune spread at all. Therefore, we postpone its de- 
tailed analysis to Sec. VI where more realistic Gaussian 
distribution will be invoked. 

A peculiarity of a single bunch beam is that its insta- 
bility is possible in a restricted region of parameters. As 
it follows from Fig. 1 at M = 1, imaginary part of the 
coefficients Vi and Vi is positive at < k — Qq < 0.5. It 
means that, without chromaticity, the instability is feasi- 
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ble only if betatron tune is located between half-integer 
and next integer (e.g. Q = 0.75 but not 0.25). This re- 
sult was first obtained in Ref. Q for a short bunch where 
the chromaticity was ignored by the model. Actually, it 
is apparent from Eq. (16) that the chromaticity is an es- 
sential factor, mainly because it triggers the head-tail in- 
teraction. For example, without chromaticity the bunch 
is quite stable at k — Q = —0.25, B = 0.5. However, the 
rigid mode becomes unstable if #oC = 1 obtaining the 
growth rate Im Av 00 ~ (0.35 — 0.26) q = 0.09 q (first 
term in this formula is the head-tail contribution, and 
second one - turn by turn interaction). Of course, chro- 
maticity of opposite sigh could prevent such a situation 
but higher modes instability would be enforced by this, 
as it follows from Fig. 6. 

The last term in Eq. (16) is a part of the collective 
interaction which describes the field variation of preced- 
ing bunches inside the considered one. Actually, only 
the nearest bunch gives a noticeable contribution to this 
part so that the result does not depend on number of 
bunches, in practice. Influence of this addition looks 
much like to the head-tail interaction which statement 
can be checked by comparison of Figs. 5 and 7. However, 
the effect strongly depends on collective mode as it is de- 
scribed by the coefficient V^k). For the example above, 
its contribution to the rigid mode instability is about 
— qoho(6oQ/16 which is less then the "usual" head-tail 
effect in order of magnitude. Probably, this part of the 
resistive wake is negligible at any conditions. 



V. LOW SYNCHROTRON FREQUENCY 

The boxcar model gives a broad outlook of the prob- 
lem but maybe it omits some important details being 
not sufficiently realistic itself. Therefore another point of 
view is developed in this section based on the approach 
Q s <C AQ which is rather characteristic of many proton 
machines. As it is shown in previous section and illus- 
trated by Fig. 3, space charge suppresses all the modes 
Yn,m<n{T, u) in this limit. Therefore the following results 
are actually concerned only with the modes Y n _ n (r) which 
will be denoted later simply as Y n (r, u). Following Ref.fl] 
with resistive wakefield added, one can show that the 
space part of this function satisfies the equation: 
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As it is shown in Ref. [f| and the eigentunes of 

Eq. (23) are ~ n(n + l)Q 2 s j AQ in order of value. There- 
fore, with rather small synchrotron frequency and for not 
very high modes, it can be simplified using the approxi- 
mation \v\ <C AQ which results in: 



U 2 (r)Y"(r)=R(r) 



(25) 




s/M^Vxj y(r') dr' - V 2 B J v{t'){t'-t) dr' 



(26) 



Boundary conditions of this equation are evident from 
the relation U 2 (±l) = which follows from definition 
(24) and will be reinforced by examples in the subsequent 
sections. Therefore, any appropriate solution of Eq. (25) 
should satisfy the relations: 



R(±l) = 0. 



(27) 



Because Eq. (25) is linear and uniform, initial conditions 
Y(l) = 1, i?(l) = can be used in practice to cal- 
culate the function R(t) everywhere with any trial u, 
and to separate thereafter the eigentunes v n assuring 
the condition R(— 1) = [H]-[6j]. The method is espe- 
cially effective at qo = to determine the basic modes. 
In particular, it confirms that Legendre polynomials are 
solutions of the boxcar bunch. Generally, it is easy to 
show that the basic solutions of any bunch are regular 
functions satisfying the orthogonality conditions: 



Y ni (T)Y n2 (T) P ( T )dT 



(28) 



Therefore, with enough small qo, additions to the eigen- 
tunes can be found by standard perturbation methods 
which way results in expression like Eq. (16): 



Av n = q a 



1\P^fn + VM(g n Vx+ih n V2B 
V #r> v 
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Eqs. (18)-(20) can be used as well with appropriate eigen- 
functions y n (i~) = Y n (r)p(r) exp(— iOoQr) to calculate 
the factors /„, g n , h n All the coefficients A„.„ = 1 in 
this case due to normalization condition (28). 



A. Gaussian bunch with low wake 

Truncated Gaussian bunch is considered in this subsec- 
tion for a comprehensive investigation of the instability. 
Its distribution function is: 
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at 
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where the normalizing coefficient C depends on a. Other 
involved functions are: 



p(r)=C[ Yexp(T 2 )erf(T)-T 



2CT 3 / 2T 2 
— 1 + — 



and 



* 2 n 



2CT £ 



2a 2 T 2 



(31) 



(32) 



where T 2 = (1 — t 2 )/(2ct 2 ), and approximate expressions 
at |t| ~ 1 are added for references. 

The case a — 1/3 (3<r truncation, C = 0.016) is ac- 
tually considered below. Six basic normalized eigcnfunc- 
tions of the bunch are shown in Fig. 8. Their eigen- 
tunes have a form v n = a n Q 2 /AQ c with the coefficients 
a n which are also presented in Fig. 8 in the brackets (for 
comparison: a n — n(n+l)/2 for the boxcar model). Here 
and further, the subindex c marks the bunch center. 

The coefficients /„, g n and g n are plotted in Figs. 9- 
12 which look much like Figs. 4-7 of the boxcar model. 
Of course, it is necessary to take into account that the 
Gaussian bunch has less rms length in comparison with 
the boxcar one (1/3 instead 1/vo), so the Gaussian 
plots should be proportionally wider. The absence of 
secondary oscillations is well explicable because of more 
smooth bunch shapes. With these reservations, one can 
assert that the boxcar model provides an adequate de- 
scription of the bunch coherent instability, at least within 
the limit of low synchrotron frequency. 
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FIG. 8: Normalized eigenfunctions of truncated Gaussian 
bunch (0th-5th modes). Corresponding eigennumbers are 
given in the parentheses. 
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FIG. 9: Functions Re/ n (0oC) of truncated Gaussian bunch. 
The right-hand parts of these even functions are shown. 
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FIG. 10: Functions Im/ n ((?oC) of truncated Gaussian bunch. 
The right-hand parts of these odd functions are shown. 



B. Expanded low wake approach 

Formally, relation (29) is applicable if the spectrum 
shifts Av n are small in comparison with distances be- 
tween the basic spectrum lines which are 

i \ Qs i i -i \ Qs 

v n +\ —v n = \<*n+i - a n)-rjr- — (n + lj- 



AQ C 



Therefore, condition of applicability of Eq. (29) is for the 
lowest mode: 



|Ai/ | < 



(33) 
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FIG. 11: Functions <?n(#oC) °f truncated Gaussian bunch. The 
right-hand parts of these even functions are shown. 



FIG. 12: Functions fo n (#oC) °f truncated Gaussian bunch. 
The right-hand parts of these odd functions are shown. 



The left-hand part of this expression is close to the insta- 
bility growth rate which is essentially less of 1 probably in 
all practical cases (e.g. | A^o| < 0.1). However, the right- 
hand part can be still less, for example 0.05 2 /0.25 = 0.01. 
The example demonstrates that a violation of condition 
(33) is quite possible occasion, especially when collective 
modes instabilities of a multi-bunch beam are examined. 
Therefore, we consider this important case in greater de- 
tail, without the assumption that the multi-bunch con- 
tribution is small. 

As a first step, we need to solve Eq. (25) at V\ — 
to know the basic modes of this case. Let us denote 
corresponding eigenfunctions and eigentunes as T„ and 
v n . Then solution of the total equation can be presented 
in the form: 

Y(t) = qoV^VM f y(T>(r') ex P HC# T') dr' 



E 

n=0 



e„T„(r) 



V — v n 



(34) 



where e n are coefficients of the expansion: 

oo 

exp(iA6» r) = e„T„(r) (35) 
It immediately results in the dispersion equation: 

oo -i 

l = q V 1 VMJ2 — / T n {T)p{T)exp{-iC0QT)dT 

(36) 

In principle, new eigenfunctions T„(r) and eigentunes 
v n could be found by the same method which was used 
for Y n (r). However, it would be a more difficult prob- 
lem to calculate the coefficients e„ because the functions 



T„(t) are not orthogonal, in contrast with Y n . Therefore 
we turn back to the approximation T = Y which is cer- 
tainly acceptable at low intra-bunch interaction and does 
not violate the overall structure of Eq. (36). It results in 
the dispersion equation 



(37) 



q VWMg n (e () 



J2 

n=o v-v n - q (y/8n/0 O fn + ih n V 2 B\/jtj 



with the same coefficients /„, g n , h n as before (Figs. 9-12 
for Gaussian bunch). With an additional condition 



go\Vi\VM< 



AQ C 



(38) 



the low wakefield approximation is totally satisfied, and 
the equation gives the same result as the earlier Eq. (29). 

Another easy case is zero chromaticity when sum (37) 
holds the only term n = because g n (0) — 6 n ,o- It 
means that solely the rigid internal mode can be excited 
without chromaticity. Because it can appear inside any 
collective mode, the resulting tune shift is: 



Aw 



9o 




/o(0) 



'MVi 



k-Qp 
M 



(39) 



This expression formally coincides with Eq. (29) at £ = 
0, but can be applied at any value of the coefficient 

Generally, series (37) contains restricted number of 
summands which conclusion follows from Figs. (6) and 
(11). In particular, one can see that the terms n = and 
1 give major contributions at |#oC| <^ 3. Eq. (37) has 
two actual solutions in this case, at least one of them be- 
ing unstable. In the extreme case when inverse of Eq. (38) 
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inequality is fulfilled, one of the eigentunes is real, and 
another tune is: 



Auj 



%Vi(k)Vm [ 5 o(0oO+Si(0oO] 



ar NR5(u) 



1 + i 



[9o(0o0 + 9i(0oO] (40) 



Rather weak dependence of this expression on chromatic- 
ity engages an attention. In relative units, the addition 
is 1 at 9 ( = and about 0.8 at 6> C = 3. However, 
it should be mentioned again that the low synchrotron 
frequency limit is considered here. Role of this factor is 
discussed in the next section. 



VI. THE INSTABILITY THRESHOLD 

It could be concluded from previous analysis that the 
boxcar is a quite adequate model for investigation of the 
bunched beam instability, and only minor and almost ob- 
vious changes are needed for more realistic distributions. 
However, it would be a premature conclusion because 
the space charge tune spread is ignored at all in the box- 
car model. Meanwhile, at certain situations the spread 
can cause Landau damping and suppress many unstable 
modes of a real bunch, as it has been shown in Ref. Q]- 
Q . The rigid intra-bunch mode is the only occasion when 
this mechanism does not work and cannot prevent insta- 
bility at any combination of parameters. 

Unfortunately, at present the problem is adequately 
covered only in the limiting cases (i> 1 and fi <C 1. In 
first case, this kind of Landau dampingreally works and 
suppresses almost all internal modes [5|. However, the 
mechanism is turned off in opposite limiting case which 
was just the subject of previous section. The only con- 
clusion can be done from these facts: this stabilization 
mechanism has a threshold character and actually begin 
to work at /i >^ 1. The goal of this section is to get more 
exact estimation of the thresholds. Because the space 
charge will be treated here as a promoting and dominat- 
ing factor, the wake field is excluded from the analysis. 

It should be reminded preliminary that Landau damp- 
ing arises when coherent frequency penetrates rather 
deeply in a region of incoherent tunes of the system. 
Then the particles which own tunes are located below 
or above the coherent tune could be exited in contra- 
phases by the coherent field, transforming its energy to 
the incoherent form (beam heating). Such a mechanism 
affects the beam decoherence and creates the instability 
threshold. 

Well known practical recommendation follows from 
this statement for coasting beams: incoherent tunc 
spread should exceed the space charge tune shift to avoid 
the instability. In principle, a wake field (e.g. the resis- 
tive wall contribution) affects this criterion, however, its 
influence is small in practice if the space charge domi- 
nates in the impedance budget. The last is just the case 
of our study. 



An additional important property of bunched beams is 
that, with a coherent frequency ui, the particles undergo 
an influence of harmonics of frequencies u> + jQ 3 where 
fl s is synchrotron frequency, and j is integer. An in- 
tense energy transfer is possible if any of these frequencies 
falls in the incoherent region. Fast look in Fig. 2 reveals 
that harmonics j — m have the most chances to do this. 
Therefore, more informative graph can be obtained from 
Fig. 2 by a transformation of each curve v n ,m{li) to the 
Vn,m — m li- The results are presented in Fig. 13 by the 
solid lines of the same color as in Fig. 2. 

Averaged over synchrotron phases incoherent tunes of 
the truncated Gaussian bunch lie in the region —I < 
u/AQ c < —0.274. Drawing corresponding boundary line 
in Fig. 13 and assuming that the coherent tunes of Gaus- 
sian bunch have about the same behavior as in the box- 
car model, one can make the conclusions: (i) the lowest 
(rigid) mode (0, 0) is unstable with any \x\ (ii) the higher 
modes (n, n) can be unstable at fj, <~ 0.5; (iii) the more 
is n the less is corresponding threshold of fj,; (iv) all 
higher radial modes like (n, m < n) are stable in any 
case. 

These conclusions are in a reasonable agreement with 
results of Ref. [5] , according which only the rigid mode 
(0,0) of Gaussian bunch is unstable at [i ^> 1. One 
can anticipate from Fig. (13) that thresholds of all other 
modes are located at \i < 0.5. It is a region where the 
low ix approximation could be fitted to refine the thresh- 
olds of Gaussian bunch. To accomplish this, we consider 
Eq. (23) with qq = but without the additional simplifi- 
cation v -C AQ C . Correspondingly, boundary conditions 
(27) should be changed by the following one: 



y'(±i) 



Qi 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

ji = Q s /AQ(0) 



FIG. 13: Transformed eigentunes: solid lines - boxcar model; 
dashed lines - Gaussian bunch with 3cr truncation. 
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TABLE I: Instability thresholds of Gaussian bunch with 
3cr truncation (lower radial modes (n,n)). 



n 





1 


2 


3 


4 


Qs/AQ C < 


CO 


0.63 


0.21 


0.13 


0.095 



Obtained eigentunes are plotted in Fig. (13) by dashed 
lines above the Gaussian incoherent boundary. The 
crossing points which are the expected thresholds of the 
head-tail modes are presented in Table 1. 

Distributions with more abrupt bunch boundaries 
demonstrate similar behavior but have higher thresholds. 
For example, three modes of parabolic bunch retain the 
stability at /i — > oo |5| . The boxcar bunch is an extreme 
case which is unaffected by this sort of damping at all. 

An important conclusion follows from this or similar 
table, concerning an influence of chromaticity in real con- 
ditions. As might be expected from Figs. 6 and 11, an 
increase of chromaticity should not affect drastically the 
instability growth rate because its main effect is a sim- 
ple switch of the bunch oscillations from a lower internal 
mode to a higher one. For example, it follows from Fig. 11 
that the most unstable internal modes of Gaussian bunch 
are: n = 1 at \0 (\ = 4, and n = 2 at 1 6* C I = 6, in both 
cases the instability growth rate being about 0.45 in used 
relative units (the head-tail contribution is neglected in 
the estimations because it is reasonably small in a multi- 
bunch beam). However, now we must take into account 
that these results were obtained at fj, = 0. Let us consider 
the case fj, = 0.4 as an another example. Then the modes 
n > 2 cannot appear being suppressed by Landau damp- 
ing. Generally, only 0th and 1st internal modes could be 
unstable in this case, and chromaticity |#oCI — 10 is su f~ 
ficient to suppress both of them that is to reach a total 
beam stability. 

Of course, Table 1 is only an estimation of the thresh- 
olds because approximate Eq. (23) lies in its foundation. 
It would be a good idea to solve more general Eq. (7) 
with arbitrary Q s and AQ C , and to use the results for 
exact instability thresholds of realistic bunches. Note 
that all eigennumbers of this equation are real. There- 
fore, the lack of regular solutions with real v at some 
combination of synchrotron frequency and space charge 
tune shift would be a sign of Landau damping. However, 
the boxcar model is the only solved case at present. 



VII. SUMMARY AND CONCLUSION 

Transverse instability of a bunched beam is studied 
with synchrotron oscillations, space charge, and wake- 
held taken into account. Resistive wall wakefield is con- 
sidered as the most universal and practically important 
case. However, many results have a common sense and, 
with small changes, can be adapted to other wakes. Box- 
car model of the bunch is extensively used in the paper 
to get a general outline of the problem in wide range of 
parameters. A realistic Gaussian distribution is invoked 
in some cases for comparison and more detailed investi- 
gations of important problems like Landau damping. 

Eigenfunctions and eigentunes of the beam are investi- 
gated with both intra-bunch and inter-bunch interactions 
taken into account. Contributions of the interactions to 
the instability growth rate are studied over a wide range 
of the parameters, including effects of the bunch length 
and chromaticity. It is shown that known head-tail and 
collective modes instabilities are the extreme cases when 
one type of the interaction certainly dominates. How- 
ever, an essential influence of the intra-bunch degrees of 
freedom on the collective instabilities is especially marked 
and investigated in detail. In particular, it is shown that 
a variability of the internal bunch modes explains why 
the instability growth rate depends on the bunch param- 
eters including space charge tune shift, synchrotron tune, 
bunch length, and chromaticity. 

It is emphasized that the space charge tune spread 
can cause Landau damping and suppress the instability 
(other than the space charge sources of the tune spread 
are not considered in the paper). The phenomenon ap- 
pears at rather large ratio of synchrotron frequency to 
the space charge tune shift, lower internal modes obtain- 
ing the stability at larger the ratio. Several modes of 
Gaussian bunch are considered in the paper, and their 
thresholds are estimated by comparison of the limiting 
cases. 
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